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We examine the effect of changing the linear polarization angle 9 of incident microwaves with re- 
spect to the dc current on radiation-induced magnetoresistance oscillations in a two-dimensional 
(2D) system within the balance-equation formulation of the photon-assisted magnetotransport 
model, considering the radiative decay as the sole damping mechanism. At an extremum the am- 
plitude of oscillatory magnetoresistance R xx exhibits a sinusoidal, up to a factor of 5, magnitude 
variation with rotating the polarization angle 9. The maximal amplitude shows up generally at a 
nonzero 9, which is dependent upon the extremum in question, the 2D electron setup, the radiation 
frequency and the magnetic field orientation. These results provide a natural explanation for the 
experimental observations by Mani et al. [Phys. Rev. B 84, 085308 (2011)], and Ramanayaka et 
al. [Phys. Rev. B 85, 205315 (2012)]. 
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I. INTRODUCTION 

Microwave induced magnetoresistance oscillation in 
high-mobility two-dimensional (2D) electron systems 
has been a subject of intensive experimental^— and 
theoreticali2r— studies over the past decades. Though a 
consensus on the period and phase of the oscillation has 
long been reached, the effect of the microwave polariza- 
tion on the amplitude of the oscillation remains one of the 
most challenging and unsolved issues since the discovery 
of microwave-induced magnetoresistance oscillations. 

Early measurement on L-shapcd Hall bars indicated 
that the period and phase of radiation-induced mag- 
netoresistance oscillations are insensitive to the relative 
orientation between the microwave polarization and the 
current^ A later experiment 13 carried out on specimens 
with a square geometry in a quasioptical setup reported a 
striking result that not only the frequency and phase but 
also the amplitude of radiation- induced resistance oscilla- 
tions and the zero resistance regions are notably immune 
to the sense of circular and linear polarizations of mi- 
crowaves. This influential result raised a big challenge 
to the existing theoretical models, in which, though no 
detailed investigation was reported, some kind of polar- 
ization dependence was believed to exist) 19 ' 22 ' 23 and thus 
expedited the emergence of different scenarios capable of 
showing polarization immunity of microwave magnetore- 
sistance rcsponsei 33 ' 34 In a recent study, Mani et aiw^ 
found a strong sensitivity in the amplitude of radiation- 
induced magnetoresistance oscillations to the relative ori- 
entation of the linear polarization with respect to the Hall 
bar axis. Particularly, more detailed measurement^ by 
rotating, by an angle 0, the polarization direction of lin- 
early polarized microwaves with respect to the long axis 
of the Hall bar electron devices, showed, at relatively low 
microwave power, a strong sinusoidal variation in the di- 
agonal resistance R xx vs 9 at the oscillatory extrema. 
And, unexpectedly, the angle for the maximal oscillatory 



R xx response under a given-power linear-polarized mi- 
crowave, which is not at 0° or 90°, appears to depend 
upon the radiation frequency, the extremum in ques- 
tion, and the magnetic field orientation. So far, there 
has been no theoretical explanation for these interesting 
observations. There is an urgent need to analyze the de- 
tailed polarization dependence from a theoretical model 
for radiation-induced magnetoresistance oscillations. 

II. MAGNETORESISTANCE UNDER 
POLARIZED RADIATION 

We deal with an isotropic 2D system of short thermal- 
ization time, consisting of N s electrons in a unit area of 
the x-y plane. These electrons, scattered by random im- 
purities and by phonons in the lattice, are subjected to a 
uniform magnetic field B — (0, 0, B) in the z direction. 
When an electromagnetic wave of angular frequency w 
illuminates perpendicularly onto the 2D plane with the 
incident electric field 

Ei(t) = E is sin(wi) + E ic cos(ujt) (1) 

at z = and a dc current flows within the plane, the elec- 
tric field inside the 2D system involves a dc component, 
Eq, and an ac component 

E(t) = E s sin(wi) + E c cos(wt). (2) 

The steady-transport state of this electron system can be 
described by the drift velocity of the electron integrative 
(the center of mass) motion, consisting of a dc part, v, 
and a stationary time-dependent part, 

v(t) = v s sin(a;£) + v c cos(cj£), (3) 

in the 2D plane, together with an average temperature 
T e , characterizing the isotropic thermal distribution of 
electrons in the reference frame moving with the center 
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of massif They satisfy the following force and energy 
balance equations^ 



N s eE + N s e{v x B) + F = 0, 
eE s F s e 

-{v s x B), 



N s mu) mui 
F r e 



eE c 

muj N„muj muj 



{v c x B), 



N s eE ■ v + S p - W = 0. 



(4) 
(5) 

(6) 
(7) 



Here, 



F = J2\ U (l\0\ 2 E 911^(0^(911, wo -n«) (8) 



is the time-averaged damping force against the electron 
drift motion due to impurity scattering, and 



S p = E 1^(911)1 



E 



nujJ^)n 2 (q\\,uJo - nuj) (9) 



is the time-averaged rate of the electron energy- 
absorption from the radiation field. In Eqs. ((8| and ((9|), 
U(qu) is the effective impurity potential, ^(qu, (2) is 
the imaginary part of the electron density-correlation 
function at temperature T e in the presence of the mag- 
netic field without the electric field, cj = q\\ • v, and 
</„(£) is the Bessel function of order n with argument 
£ = [(qfii • v s ) 2 + (<j»|| -v c ) 2 ]2/uj. Note that, although con- 
tributions of phonon scattering to F and S p are neglected 
in comparison with those of impurity scattering at the 
considered low lattice temperature, it provides the main 
channel for electron energy dissipation to the lattice with 
a time-averaged energy-loss rate W, having an expression 
as given in Ref.l22l. 

The ac components v s and v c of electron drift veloc- 
ity should be determined selfconsistently in terms of the 
incident ac field Ei(t) by the electrodynamic equations 
connecting both sides of the 2D system, taking into ac- 
count the scattering- related damping forces F s and F c *£ 
However, for high-mobility systems at low temperatures, 
the effects of these scattering-related damping forces are 
much weaker in comparison to those of radiative decay^i 
and thus negligible, whence v s and v c are in fact directly 
given from Eqs. ([5]) and (J6)) by the high-frequency electric 
field E(t) inside the 2D electron system. On the other 
hand, by solving the Maxwell equations connecting both 
sides of the 2D electron gas which is carrying the sheet 
current density N s ev(t), E(t) is determined by the in- 
cident fields E[ S and E\ c based on the setup of the 2D 
system in the sample substrate^ 

If the 2D electron gas locates within a thin layer un- 
der the surface plane at z = of a thick (treated as 
semi-infinite) semiconductor substrate having a refrac- 
tive index, n s , the ac field E(t) driving the 2D electrons, 
which equals the sum of the incident and the reflected 



fields at z — and equals the transmitted field (the field 
just passes through the 2D layer), can be expressed as^i 



E(t) 



N s ev(t) 



2n 



(n + n s )eoc n + n 



(10) 



Here no is the refractive index of the air and c and eo are, 
respectively, the light speed and the dielectric constant 
in vacuum. If the 2D electron gas is contained in a thin 
layer suspended in vacuum at the plane z — 0, then 



2e c 



(11) 



By combining Eq. (fTUj) or Eq. (fTI) with balance Eqs. (0 
and (|6]), the ac drift- velocity components v s and v c , and 
thus the argument £, can be determined by the incident 
field components E- ls and E lc . 

The imaginary part of 2D electron density correlation 
function in a magnetic field, II2 {q\\ , /?), can be written in 
the Landau representation as^ 

tfafeu, O) = ^fT E C n , nl {l 2 B q\/2)n 2 {n, n', O), (12) 

n,n' 

n 2 (n,n',0)=-~ Jde [f(e) - f(e + Q)] 

x ImG„(e + Q)lmG n >{e), (13) 



where I b = y/T/\eB\ is the magnetic length, C n ,n+i(Y) = 
nl[(n + l)\]- 1 Y l e- Y [L l n (Y)} 2 with L l n (Y) being the asso- 
ciate Laguerre polynomial, /(e) = {exp[(e— ^)/T e ] + l} _1 
is the Fermi function at electron temperature T e , and 
ImG n (e) is the density-of-states (DOS) function of the 
broadened Landau level n. 

The Landau-level broadening, which results from im- 
purity, phonon, and electron-electron scatterings, is as- 



sumed to have a Gaussian form [e n = (n 
center of the nth Landau level, n = 0, 1, 2, . 



is the 



ImG„(e) = -(27r) = r 1 exp[-2(e - e n f / T 2 }, (14) 

with a B i -dependent half width expressed as 

r = (2 Wc /vrr s )5, (15) 

where r s , the single-particle lifetime or quantum scat- 
tering time in the zero magnetic field, is assumed to be 
related to the transport relaxation time Tt r or the zero- 
field linear mobility /io using an empirical parameter a 
by 22 



1/t s = 4a/r tr = 4ae/m/i - 



(16) 



For an isotropic system where the frictional force F is 
in the opposite direction of the drift velocity v and the 
magnitudes of both the frictional force and the energy- 
dissipation rate depend only on v = \v\, we can write 
F(v) = F{v)v/v and W(v) = W(v). In the Hall con- 
figuration with velocity v in the x direction v = (v, 0, 0) 
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FIG. 1: (Color online) Magnetoresistivity R xx , electron tem- 
perature T e and energy absorption S p are plotted versus 
uj/uj c under the irradiation of linearly polarized (LP) mi- 
crowaves of frequency 40-GHz and incident power Pi u — 
17.5 W/m 2 at several different polarization directions (6 — 
0°, 30°, 60°, 90°, 120° and 150°) for a system described in the 
text at temperature T = 1.5 K. 
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FIG. 2: (Color online) Magnetoresistivity R xx induced by 
linearly polarized microwave of frequency 40 GHz and incident 
power = 17.5 W/m 2 at 6 = 30° is shown for positive and 
reverse magnetic fields over the range —2.5 < uj c /uj < 2.5 (a). 
The 9 dependence of R xx at extrema P + l, V + l, P + 2, and 
V+2 (b) and at extrema P~l, V~l, P~2, and V~2 (c). 



or the current density J x — J — N s ev, and J y = 0, the 
longitudinal linear (v — > 0) resistivity under the incident 
radiation of Eq. ([T]) due to impurity scattering, can be 
written in the form 

^ = -E9^|^ E ^(Offafol.nw), (17) 

qy s n— — oo 

where IlJ,{q h 0) = dn 2 (q\\, Q)/dQ. 



III. OSCILLATORY RESISTIVITY VERSUS 
POLARIZATION DIRECTION 

We consider an incident radiation of form (JTJ linearly 
polarized along the direction having an angle with 
respect to the x-axis: E- ls — (Ej U cos 6, E- 1U] sin 6) and 
Ei C = 0. The 2D electron gas, having carrier sheet den- 
sity N s = 2.2 x 10 15 m -2 and zero-temperature linear mo- 
bility /in = 800m 2 /Vs from short-range impurity scat- 
tering in the absence of the magnetic field, is assumed 
to locate within a thin layer under the surface plane at 
z = of a thick GaAs-based substrate with a refractive 
index of n s — 3.59. 

Figure 1 shows the calculated energy absorption rate 
Sp, the electron temperature T e and the longitudinal lin- 
ear resistivity R xx of this system irradiated by linearly 
polarized (LP) incident microwaves having frequency 



lj/2tt = 40 GHz and amplitude = 1.15 V/cm (i.e., in- 
cident power P iu = 17.5 W/m 2 ) at different polarization 
directions (6 = 0°, 30°, 60°, 90°, 120°, and 150°), together 
with the dark resistivity, as functions of the normalized 
inverse magnetic field uj/lo c (uj c — eB/m). The lattice 
temperature is assumed to be T = 1.5 K and the Landau 
level broadening parameter is taken to be a = 3. 

We see that the electron energy absorption S p and 
thus the electron temperature T e , both showing a marked 
main peak at cyclotron resonance lj/uj c — 1 and sec- 
ondary peaks around its harmonics uj /lo c = 2, 3, 4, are 
essentially the same for all different polarization direc- 
tions of radiation. This can be understood directly from 
expression §§§ of S p in the case of relatively low strength 
of incident microwave when the dominant contribution 
to it comes from the terms n = ±1 and J 2 _ 1 (^) ~ £ 2 . 
Writing explicitly the ^-dependent expression of £ 2 in the 
case of incident plane-polarized radiation having an angle 
9 with respect to the x-axis we can see that, after sum- 
ming over all the directions of t/|| , the angle dependence 
disappears due to sin 2 8 + cos 2 8 = 1. Furthermore, the 
energy-loss rate W has ^-dependent behavior similar to 
that S p . Thus the electron temperature T ei determined 
by the energy-balance equation (7) involving only S p and 
W, has similar ^-dependent behavior. The situation is 
different for the frictional force F and the resistivity R xx , 
because of additional qy-direction dependent weighted 
factors showing up inside the qn summation in the ex- 
pressions ((5|) and (JT7J) , leading to sensitive ^-dependence 
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of R xx as seen from the numerical results shown in the 
upper part of Fig. 1. The linear-polarized microwave ex- 
cited magnetoresistivity R xx oscillates strongly, having 
the same oscillatory period and nearly the same phase 
and, at cyclotron resonance uj/uj c = 1 or its harmonics 
lo/u c = 2,3,4,..., taking equal values for all different po- 
larization directions 9. We can label, for all 9, the main 
peaks and valleys of the oscillatory resistivity respectively 
as PI, P2, and VI, V2, etc. The amplitude of resistiv- 
ity oscillation, however, varies sensitively with changing 
the polarization direction of the microwave. The max- 
imal or minimal amplitude of a peak or a valley shows 
up at different polarization angles 9 for different peaks 
or valleys. For instance, the maximal and minimal am- 
plitudes show up, respectively, at 9 w 59° and 9 « 149° 
for PI, at 9 « 25° and 9 w 115° for VI, at 9 « 12° and 
9 w 102° for P2, at 9 « 8° and 9 w 98° for V2, and show 
up at polarization angles closer to 9 = 0° and 9 = 90° 
for higher order peaks and valleys. The exact polariza- 
tion angle for the maximal or minimal amplitude of an 
cxtrcmum appears to depend on the radiation frequency 
and the 2D electron setup in the substrate. In obtaining 
the above results the sample setup is so assumed that 
2D electrons locate within a thin layer under the surface 
plane at z = of a thick semiconductor substrate having 
a refractive index of n s = 3.59, and the relevant damping 
in question, i.e., the radiative decay, is fully determined 
by this sample setup from the electrodynamic equation 
(|10p . Except for this, no other damping parameter nor 
any mechanism capable of producing polarization rota- 
tion is introduced in the present model. 

Next we present the results obtained under magnetic 
field reversal. Figure 2(a) plots the longitudinal resistiv- 
ity R xx vs u) c /uj over the range —2.5 < lj c /uj < 2.5 for the 
system irradiated by a linearly polarized microwave along 
the 9 = 30° direction with frequency gj/2it = 40 GHz and 
incident power p w = 17.5 W/m 2 . The extrema of inter- 
est here are labeled as P + l, V + l, P + 2, and V + 2 for 
those in the domain B > 0, and P~l, V~l, P~2, and 
V~2 for those in the domain B < 0. The values of lon- 
gitudinal resistivity R xx at P+l, V+l, P+2, and V+2 
are plotted in Fig. 2(b), and those at P~l, V~l, P~2, 



and V~2 qre plotted in Fig. 2(c), as functions of the 
microwave polarization direction 9. The present treat- 
ment is for an isotropic system. By symmetry we al- 
ways have P+ 1(9) = P+l (tt + 6) , P~ 1 (6) = P~ 1 (tt + 6) , 
P+l(9) = P~1(tc — 9), etc. If the maximal amplitude of 
an extremum, e.g., P+l, shows up at 9 ^ 0, the maximal 
amplitude of the corresponding extremum at the reverse 
magnetic field, P~l, must be at a different polarization 
angle. The effect of asymmetry in a real sample itself 
would produce further complexity. 



IV. SUMMARY 

We have examined the effect of changing the polariza- 
tion angle 9 of the incident linearly-polarized microwaves 
with respect to the dc current on radiation-induced mag- 
netoresistance oscillations within the balance-equation 
formulation of the photon-assisted magnetotransport 
model. 

The present investigation takes the incident microwave 
field, rather than the ac field inside the 2D electron sys- 
tem, as the input quantity, allowing a direct determina- 
tion of the dominant damping mechanism in the high- 
mobility system, the radiative decay, from the experi- 
mental sample setup without introducing any artificial 
damping parameter, thus enabling a direct comparison 
with experimental measurement. 

We find that the amplitude of radiation-induced mag- 
netoresistance oscillation varies sensitively with changing 
9. At an extremum the amplitude of oscillatory magne- 
toresistance R xx exhibits a sinusoidal, up to a factor of 
5, magnitude variation with rotating the polarization an- 
gle 9. The maximal amplitude shows up generally at a 
nonzero 9, which is dependent upon the extremum in 
question, the 2D electron setup, the radiation frequency 
and the magnetic field orientation. 
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